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O ' Abstract 
(N 

^ ■ We consider a cooperative wireless network in the presence of one of more eavesdroppers, and 

exploit node cooperation for achieving physical (PHY) layer based security. Two different cooperation 
schemes are considered. In the first scheme, cooperating nodes retransmit a weighted version of the source 
. signal in a decode-and-forward (DF) fashion. In the second scheme, while the source is transmitting, 

cooperating nodes transmit weighted noise to confound the eavesdropper (cooperative jamming (CJ)). 
We investigate two objectives, i.e., maximization of achievable secrecy rate subject to a total power 
constraint, and minimization of total power transmit power under a secrecy rate constraint. For the first 
design objective with a single eavesdropper we obtain expressions for optimal weights under the DF 



C/2 



> 

OO . protocol in closed form, and give an algorithm that converges to the optimal solution for the CJ scheme; 

m ■ 

while for multiple eavesdroppers we give an algorithm for the solution using the DF protocol that is 
_ guaranteed to converge to the optimal solution for two eavesdroppers. For the second design objective, 

, existing works introduced additional constraints in order to reduce the degree of difficulty, thus resulting 



X 



in suboptimal solutions. In this work, either a closed form solution is obtained, or algorithms to search 
for the solution are proposed. Numerical results are presented to illustrate the proposed schemes and 



H ' demonstrate the advantages of cooperation as compared to direct transmission. 

Index Terms 

Secrecy rate, node cooperation, physical layer based security, semi-definite programming. 
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I. Introduction 

Privacy and security issues play an important role in wireless networks. Although security is typically 
addressed via cryptographic approaches, there have been several attempts at addressing security at the 
physical layer, following the pioneering work of |8|. Recently wireless physical (PHY) layer based 
security from a information-theoretic point of view has received considerable attention, e.g., 131-1171. The 
wiretap channel, first introduced and studied by Wyner m, is the most basic physical layer model that 
captures the problem of communication security. Wyner showed that when an eavesdropper's channel 
is a degraded version of the main channel, the source and destination can achieve a positive perfect 
information rate (secrecy rate). The maximal secrecy rate from the source to the destination is defined as 
the secrecy capacity and for the degraded wiretap channel is given as the difference between the rate at the 
legitimate receiver and the rate at the eavesdropper. The Gaussian wiretap channel, in which the outputs at 
the legitimate receiver and at the eavesdropper are corrupted by additive white Gaussian noise (AWGN), 
was studied in [9|. Along the same lines, the Gaussian MIMO wiretap channel was investigated and the 
secrecy capacity of the MIMO wiretap channel was established in terms of an optimization problem over 
all possible input covariance matrices ifTOl . |[T3l . There have also been some recent works focusing on 
secrecy rates based on partial CSI or channel statistics 131, HI, 1151 . In |[T5l . the authors derived the 
ergodic secrecy capacity of Gaussian MIMO wiretap channel and showed that a circularly symmetric 
Gaussian input is optimal. 

The secrecy rate is affected by channel conditions between the source and the destination and also 
channel conditions between the source and the eavesdroppers. A low cost approach to increase the 
achievable secrecy rate by exploiting/mitigating channel effects is node cooperation via relays |[T6l - ll2ri . 
A two-stage cooperative approach was recently proposed in l22ll . ||231 . and their extended version l24l . 
In l24l . the source first transmits locally to a set of trusted relays, and subsequently, the relays retransmit 
a weighted version of the signal that they heard (amplify-and-forward (AF)), or a weighted version of 
the decoded signal (decode-and-forward (DF)). Alternatively, the relays can transmit weighted noise to 
confound the eavesdropper while the source is transmitting (cooperative jamming (CJ)). In all cases, the 
objective is to select the weights so as to maximize the secrecy rate under total power constraints, or to 
minimize the total power under a secrecy rate constraint. The results in l22l - l24ll contain sub-optimal 
weights for both a single eavesdropper and multiple eavesdroppers, due to the difficulty of solving 
the associated optimization problems. In particular, several criteria for sub-optimal weight design were 
proposed such as completely nulling out the message signal at all eavesdroppers for DF and AF, and 
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completely nulling out the jamming signal at the destination for CJ. These sub-optimal weights may yield 
a reduction in the achievable secrecy rate or minimization of total power. 

In this paper, we consider the same scenario and problem as in |[22l-| |24l . but focus on obtaining the 
optimal solution for the DF and CJ schemes. Obtaining the solution for the AF scheme is a more difficult 
problem and will be addressed in future work. Exploiting certain properties of the objective functions 
and the constraints, enables us to either obtain closed form solutions, or, if a closed form solution is not 
possible, propose algorithms to search for the solution. 

The remainder of this paper is organized as follows. The mathematical model is introduced in ^ 
In we derive the optimal relay weights that maximize the achievable secrecy rate subject to a total 
power constraint in the presence of a single eavesdropper for the DF and CJ protocols, and multiple 
eavesdroppers for DF protocol. In ^TV] we study the optimal weights that minimize the total power under 
a secrecy rate constraint for the DF and CJ protocols. Numerical results in ^ are presented to illustrate 
the proposed solutions. Finally, WII provides some concluding remarks. 

A. Discussion of related work 

Our work falls under the general scenario where the source-destination communication is aided by 
a relay or a helper. Relevant results include the work of lfl6l . where multiple users communicate with 
a common receiver in the presence of an eavesdropper, and the transmit power allocation policy is 
determined that maximizes the secrecy sum-rate. In |17|, a source, destination, eavesdropper and relay 
model is considered, in which the relay transmits a noise signal in order to jam the eavesdropper. The 
rate-equivocation region is derived to show gains and applicable scenarios for cooperation, with the 
equivocation denoting the uncertainty of the eavesdropper about the source message. 

A generalization of I1I61I and ifTTl was proposed in |[T8l . in which the helper transmits signals from 
another source encoder. In fT9l, inner and outer bounds on the rate-equivocation region were derived for 
the four-node model for both discrete memoryless and Gaussian channels. In lIlTl . the secrecy rate of 
orthogonal relay and eavesdropper channels was studied. 

Our work in this paper is different from the aforementioned works in the sense that we address the 
more general case of multiple relays and multiple eavesdroppers. Also, existing works primarily focus 
on rate-achieving relaying strategies. In our work, we consider pre-defined cooperative schemes without 
claiming that those schemes are optimal, and determine relay weight and power allocation design that 
optimize the achievable secrecy rate subject to a power constraint, or minimize the transmit power subject 
to a secrecy rate constraint. 
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B. Notation 

Upper case and lower case bold symbols denote matrices and vectors, respectively. Superscripts *, T 
and t denote respectively conjugate, transposition and conjugate transposition. Tr(A) denotes the trace 
of matrix A. A ^ means that A is a Hermitian positive semi-definite matrix. rank(A) denotes the 
rank of matrix A. ||a|| denotes Euclidean norm of vector a. I.„ denotes the identity matrix of order n (the 
subscript is dropped when the dimension is obvious). f'{x) and f"{x) denote first- and second- order 
derivatives of f{x), respectively. 



Source 




Fig. 1. System model: source S wishes to communicate to destination D in the presence of J eavesdroppers, Ei,. . . ,Ej. 
The relays, Ri, . . . , Rm implement decode and forward (DF) or cooperative jamming (CJ), based protocols. In each case 
the objective is to select the relay weights (w) and the source power (Ps) to maximize the achievable secrecy rate subject to a 
total power constraint (Po), or to minimize the total power constraint under a secrecy rate constraint. 

II. System Model and Problem Statement 

We consider a wireless network model depicted in Fig. [T] consisting of one source node S, a set 
of relay nodes {Ri, i = 1,. . . ,N), a destination node D, and a set of J passive eavesdroppers 
{Ej, j = 1, . . . , J). The symbols used in the paper are listed in Table U 

The source message is uniformly distributed over the message set W = {1, 2, 2"^}, which is 
transmitted in n channel uses. Here, R denotes the source rate (unit: bits per channel use) and the 
message has entropy nR bits. A stochastic encoder at the source maps each message to a codeword 
from an alphabet of length-n. For the purpose of evaluating the achievable secrecy rate, we assume that 
the codewords used at the source are Gaussian. We consider a time division multiple access system, in 
which there are n time units in each transmission slot. In a time unit, the average power of an encoded 
source symbol is normalized to unity. The noise at any node is assumed to be zero-mean white complex 
Gaussian with variance o"^. Each node is equipped with a single omni-directional antenna and operates 
in a half-duplex mode. 
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TABLE I 

Mathematical notation 



N 


number of relays 


J 


number of eavesdroppers 


Po 


total power (source power plus the relays' power) 


Ps 


transmit power at the source 




noise variance 


ho 


baseband complex charmel gain between the source and the destination 


hi 


baseband complex channel gain between the ith relay and the destination 


tti 


baseband complex channel gain between the source and the ith relay 


aoj 


baseband complex channel gain between the source and the jth eavesdropper 


Qij 


baseband complex channel sain between the ith relay and the ith eavesdropper 


W 


weight vector at the relays, [wi, ■ ■ ■ , tiijv]"^ 


h 


[fti, ■ ■ ■ , KnY' 








hht 


R^, 




Notes: 


The index j is dropped when J = 1. 



All channels are assumed to be flat fading. We assume that global channel state information (CSI) is 
available, including the eavesdroppers' channels. This corresponds to the cases where the eavesdroppers 
are active in the network and their transmissions can be monitored |[25l . We should note that there have 
been some recent works focusing on secrecy rates based on partial CSI or channel statistics (e.g., chapter 
5 in [7], and [4|). Adapting the proposed work to cooperative schemes that uses partial CSI or channel 
statistics will be considered in future work. 

Similarly as in |[20l . |[24ll . we assume that the source encoding scheme, the decoding methods at 
destination and eavesdroppers, and the cooperative protocol, are all public information. 

Let us fix the relay weight vector w and the source transmission power Ps. Then the following 
expressions give the rates with the destination and the eavesdroppers as a function of w, Pg, the noise 
cr^, and the various channel gains. For the DF-based protocol, the rate at the destination and the jth 
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eavesdropper are, respectively 

i?d = -log(^l + J, (1) 

Ri = ^logh+ (2) 

where the scalar factor 1/2 is due to the fact that two time units are required in two stages. Here we have 
assumed an additional constraint, i.e., Pg > P™™ where p™™ is the minimum source power requirement 
for cooperative nodes to correctly decode the source message with high probability. We assume p™'" 
is known a priori. For the CJ-based protocol, the rate at the destination and the jth eavesdropper are, 
respectively 




(3) 



(4) 



wtR^w + a"^ J 

For both the DF and CJ protocols, the achievable secrecy rate in the presence of J eavesdroppers is given 
by 1261 

Rg = max{0, — max Rl}. (5) 
In particular, when J = 1, i.e., a single eavesdropper, the secrecy rate in ^ becomes llT3l . |[T4l 

Rs = max{0,Rd- Re}. (6) 

We consider the practical case in which the system can be designed so that the secrecy rate is positive. 
In that case, the achievable secrecy rate can be rewritten as 

Rs = Rri— max Ri. (7) 
i<i<J 

Achievability of the above rates can be shown based on existing results for MIMO wire-tap channels, 
such as Ifm . |[T4l . |[T2 I. for one eavesdropper, and 11261 for multiple eavesdroppers. 

The CJ scheme can be viewed as a 1 x 2 SIMO system, so that MIMO results are directly applicable. As 
discussed in |[24ll . for the DF scheme, MIMO results are applicable if we assume that the received signal 
at destination/eavesdropper at time i depends only on the relays' transmitted encoded signals at time i 
(though a relay's transmitted signal at time i depends on its received signal before time i). This is usually 
referred to as the "memoryless relay channel" ifTTl . |[T9l . For convenience, as in ll24l we focus on case in 
which the relays use the same codewords as the source to re-encode the signal before transmission; in that 
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case the rates of source-destination and source-eavesdropper links admit simple closed-form expressions 

lEl. 

The problems addressed in this paper are described as follows. 

Problem 1 (maximize secrecy rate under power constraint using DF): Given total power (source 
plus relays) Pq, select source power Pg and cooperative nodes' weights w to maximize the secrecy rate 

max Rs s.t. + ||w||2 = Pq, P, > P™". (8) 

where Rg is given by ([T]), ([2]) and ([T]). 

The solution of Problem 1 for one eavesdropper is provided in Section IIII-A1[ and for multiple eaves- 
droppers in Section Ull-B I 

Problem 2 (maximize secrecy rate under power constraint using CJ): Given total power (source 
plus relays) Pq, select source power Pg and cooperative nodes' weights w to maximize the secrecy rate 

max Rg s.t. Pg + llwlP = Pq. (9) 

where Rg is given by (l3]l, dH) and dV). 

The solution of Problem 2 for one eavesdropper is provided in Section IIII-A2I 

Problem 3 (minimize transmit power under secrecy rate constraint using DF): Given the secrecy 
rate constraint R^, select the source power Pg and cooperative nodes' weights, w, to minimize the total 
power (source plus relays) 

min [Po = Ps + ||w||2] s.t. Rg = Pg > P^^ (10) 

where Rg is given by ©, Q and 

The solution of Problem 3 is provided in Section IIV-AI 

Problem 4 (minimize transmit power under secrecy rate constraint using CJ): Given the secrecy 
rate constraint R^, select the source power Pg and cooperative nodes' weights, w, to minimize the total 
power (source plus relays) 

mm \Po = Pg + \\y^f] s.t. P, = P^ (11) 

where Rg is given by (l3]l, dH) and 

The solution of Problem 4 is provided in Section IIV-BI 

Before proceeding, we provide Lemma [U and |2l which will be the basis of the results to follow. Please 
see Appendix |A] and |B] for details. 
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Lemma 1: Let r and s be (known) linearly uncorrelated vectors. Let 6 be the argument o/r^s. Denote 
i = \/— T. The matrix rr^ — ss^ has only two nonzero eigenvalues, i.e., r/i > and r]2 < 0, given by 

rji = ||r||^ — |c2||r^s|, r/2 = ||r||^ — |c4||r'''s|. (12) 

The corresponding eigenvectors are 

ei = ci(r + |c2|e^("-^)s), eg = C3(r + |c4|e^("-^)s) (13) 

where ci = + |c2P||sP — 2|c2||rts|, C3 = l/-\/||rP + |c4p||sP — 2|c4||rts| , 2|c2||r^s| = 

IklP + - -v/dkP + l|sP)^ - 4|rts|2, and 2|c4||r'l's| = ||rp + ||sp + -Y/dlrp + ||s||2)2 - 4|rts|2. 

Lemma 2: Let di and d2 be (known) unit-norm vectors. Let (j) G (— 7r,7r] be the argument 0/ d|di, 
r = |djd2|, and consider < q < 1. The solution of 

max zM2d2Z s.t. zMidtz = q, ||z|| = 1 (14) 

z 

/5 g/ve?i Z?}' z° = cidi +C2d2, where C2 = ^y{^^^'q)/lY^^r^ and ci = (rc2 — ^)e^^'^~'f'\ The maximum 
is z°M2d5z° = 1 - (r^r^- y/{l-r^)qf. 

III. Secrecy Rate Maximization under Power Constraint 
In this section, we address Problems [1] and |2l 



A. One Eavesdropper (J = 1) 

In this subsection we address Problem [T] for the case of a single eavesdropper. Since J = 1 we drop 
the superscript j from Ri. 

1 ) DF-based protocol: Problem [T] can be recast as 

1 /a2 + /y^ioP + wtR^\ 

max -log — — 77; r:- (15) 



s.t. R>Pn™",wtw = Pn-P,. 



-(j--, WW = /-o 

Denote ui = g/||g||, U2 = h/||h|| and ( = |u|u2|. Let 9 G (— 7r,7r] be the argument of UgUi. The 
solution of Problem [T] is given in the following theorem. The proof is given in Appendix O 
Theorem 1: The solution of di5D is given by 



po^. pt"" if JiPt"") > jm, ^^^^ 

Pn else; 
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„ C1U1+C2U2 if J(Po"^-) > J(Po); 

w = < (17) 
else 



where the function J{Ps) is defined by 



+ fil/iQp + (Pq - fi)||hf L(z(P,)) 



and 



L{z) ^ 1 - {CVT^z - ^/{l-C,^)zf (19) 
z(P,) = {B- - 4C)/(26C) - a/h (20) 

a = + PsboP, 
6 = (Po-P.)||gf, 

c = a^ + P,|/io|2 + (Po - Ps)||hf(l - (2), 
d=(Po-P.)||hf(l-2C'), 



/ = 2(Po-P,)||hfCv/l-C2, 

B = {2a + h)/{a^ + ah), 

f{2a + 6)2 + 4(a(i + 6c)2 



C 



^P{o? + a6)2 + 4(a(i + 6c)2(a2 + a6) 



C2 = y^W^ V(l - ^(n°))/(l - C^) 

ci = V^o - P° (Cc2 - .MP^)e<^-'\ 

Remarks: In particular, Theorem [T] states that, depending on the relative values of J(P™°) and J(Po), 
the optimal solution is either i) the source uses minimum power and the relay weights are a linear 
combination of the normalized relay and eavesdropper channel vectors, or ii) the source uses all the 
power and the relays are unused. 

2) CJ -based protocol: Problem |2] is recast as 

max log f 1 + ^^^) - log (l + -F^^) (21) 

s.t. wtw = Po-P«,Ps G [0,Po]. 

By denoting = \/Pq — Pg x, vi = h/||h||, V2 = g/||g||, the problem of ((2T]) can be rewritten as 

, A ^ Ps\ho? ^ , / ^ Ps\go? \ 

max log 1 H 1 - log 1 H : (22) 

V (Po-P.)||hPxtvivIx + a2y (Po-P,)||gPxtv2vTx + a2y 

s.t. xtx = 1,P, G [0,Po]. 
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Let X be a feasible point. Denote xWiv|x = z, z ^ [0, 1]. For fixed z, a larger XW2V2X results in a 
larger objective value. With this and from Lemma |2l we know that the optimal XW2V2X equals G{z) 
where G{z) = 1 — {rj^/l — z — — rf)z)'^, rj = |v|v2|. With these, we can rewrite the optimization 
of ^ as 

max log f 1 + 7- ) - log f 1 + 7- , ) (23) 

Ps,z ""V {Po - Ps)aiz + a2 J iPo-Ps)a3G{z) + aJ 

s.t. [o,i],P, e [o,Po] 

where ai = ||h||V|/ioP, "2 = (T^/\ho\'^, 03 = ||g||VboP and 04 = aVlsoP- 

The problem of ( |23l ) makes sense when the maximum is greater than zero, i.e., when positive secrecy 
rate is achieved. The conditions under which positive secrecy rate is achieved are given in the following 
lemma. The proof is given in Appendix |Dl 

Lemma 3: The condition under which positive secrecy rate is achieved is: 

• a2 < 04 (i.e., |/iop > 1^0 P)/ 

• a2 > 04, -Po > ("2 — 04) / i(^3G{zo) — aizo) where zq is the unique root of a^G'^z) = ai given by 
zo = 1/(1 + ul), uo = [ai/as - (2r?2 _ l) + ^(ai/aa)^ + 1 - 2(ai/a3)(2r?2 _ l) ]/{2ij^l - r/2). 

Remarks: The second condition means that if the source-destination channel is weaker than source- 
eavesdropper channel, direct transmission can not achieve positive secrecy rate, at that time, relays should 
be used and the total power (source plus relay) should be greater than a threshold. 

In the following analysis, we assume the conditions in Lemma [3] hold. Denote the objective in (|23] ) 
as Mi{z). It is easy to show that M[{z) > 0. Thus, z = is not the optimal point. Before proceeding, 
we give a suboptimal solution which turns out to be the suboptimal solution proposed in [i24il which is 
a special case corresponding to z = (please see Appendix |E] for details). 

Lemma 4: When z = is fixed, a suboptimal solution is given: 

. ifa2> 04, Po < ("2 - o;4)/(a3(l - r/^)), then P^.sub = 0; 

• '/ {"2 < "4 or a2 > 04, Po > ("2 - a4)/(a3(l - v"^))} (Pq + 04)04 > (Pq + a2)(Poa3(l - 
rj^) + 04), then P^^^ub = -Pq." 

• if {"2 < "4 or a2 > 04, Po > ("2 - a4)/(a3(l - V"^))} (Pq + 04)04 < (Pq + 02)(Po03(l - 

■rf) + 04), then Pg^snh = (-£3^2 + \/ cldl - 03^2(1 - ci2)(a2 - C3))/(d2(l - ^2)), where d2 = 

03(1 - -rf), C3 = 04 + Po03(l - -rf). 
Now we proceed. The methodology to solve the problem of (1231 ) is: 1) fix z, find the optimal P^; 2) 
fix Ps, find the optimal z. Based on this, we propose an algorithm to search for the optimal Pg and z as 
follows. 
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Algorithm 1: Take a feasible point z^^^ as initial point. Subsequently, find the optimal pj^^ and then 
the optimal z''^\ Then find the optimal Pp\ and so on. The procedure converges to the optimal P° and 

z°. 

The algorithm [T] is not complete without providing the methods to find the optimal z for fixed Ps 
and the optimal P3 for fixed z. Next, we provide such methods. First, we consider the problem: find 
the optimal Pg for fixed z. This corresponds to an optimization problem of a single variable Pg, and 
the maximum is achieved at either 0, Pq, or the points with zero derivative. Setting the derivative of the 
objective to zero leads to the following quadratic equation 

AiP2 + BiP, + = (24) 

where Ai = aidi{l — di) — 6iCi(l — 61), Bi = 2aici(di — 61), Ci = aici(ai — ci), oi = PqOiz + 02, 
bi = a\z, c\ = PoasG{z) + a4 and di = asG{z). We can obtain Ps explicitly as a function of z. Second , 
we consider the problem: when Ps is fixed, find the optimal z. This corresponds to an optimization 
problem of a single variable z, and the maximum is achieved at one of the following points: 0, 1 and 
the points with zero derivative. The problem to solve can now be rewritten as 

max Rs{z) = log ■ log — (25) 

z a\z + 03 OL'iGyz) + C3 

s.t. z G [0, 1] 

where 03 = a2/(Po - Ps), 63 = ("2 + Ps)/(Po - Ps), C3 = 04/(^0 - Ps) and ^3 = (04 + Ps)I{Pq - Ps). 
The derivative of Rs{z) given by 

, / s A ^ ^ -P^ ( azG\z) oi \ 

dz Po-P, \^(a3G(z) + C3)(a3G'(z)+d3) {a^z a^){a^z ^ h^) ) ' 

It is easy to verify that ^^(0) > 0' ^s(l) < 0- Thus, the optimal z must be the points with zero 

derivative. Note that Rs{z) > holds only when aiz + < a3G{z) + C3, which determines an interval 

{z,z) C [0, 1]. Here z and z can be expressed in closed form from the fact: if G{z) = fiiz + has 



real roots over [0,1], then its roots can be expressed as z = 1/(1 + Uq) where uq = (r/^1 — rp' ± 



^J^rf■{l — ff) + (/32 — 1 + 'rf'){'rf' — Pi — 132)) / (/32 — 1 + ??^). With these, we can restrict our attention to 
the root of ( [261) over {z, z). 

To proceed, we need the following result. The proof is given in Appendix [0 

Property 1: For Rs{z) defined in (l25l ). ^4^1 , < Ofar the stationary point z' G {z,z) (i.e., the point 



dz^ 

I — n\ 

dz 



with ^1 , = Oj. 
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According to Property \T\ we know that the equation (l26l ) has a unique root z' such that when z < z', 
> and when z > z', < 0. This property ensures that the Newton method would be very 
effective in searching for z' and would enjoy quadratic convergence. 



B. Multiple Eavesdroppers 

We now turn to the case of multiple eavesdroppers (J > 1). We restrict our attention to the DF protocol. 
The CJ protocol case with multiple eavesdroppers is a more difficult problem and will be addressed in 
future work. 

Problem [T] now becomes 

. 1 a2 + P,|/ioP + wtR;,w 

max mm — ioe; ■. — (27) 

P.,w jei 2 ^ + P,\gQj\-2 + wWgW 

s.t. Ps G [Po"'", Po], wtw = Po-Ps 

where / = {I,-- - , J}, is defined as in the case of a single eavesdropper. By letting w = 



-y/^j^-^x, we can rewrite the above problem as 



max mm — iog ■. — (28) 

p.,x jei 2 cT2 + P,|5o,f + (Po-P.)xtR^x 

s.t. G [Po™",Po], xtx = l. 

Before proceeding, we give a suboptimal solution which turns out to be the suboptimal solution in 124]. 
In |[24l . if > J+ 1, a suboptimal solution is obtained when an additional constraint is added: w^^G = 
where G = [gi,g2,-'" )gj]- This additional constraint means nulling the energy to the eavesdroppers, 
and this requires N > J so that we have enough degrees of freedom to ensure this is possible. The proof 
is given in Appendix [Gl 

Lemma 5: When the constraint w^G = is added, a suboptimal solution is obtained as 

_ Po if /2(Po) > /2(Po-^); 







where the function /2(Ps) is defined by 



/2(,Psj - 2 I P T\ 121 ^^^> 

+ P^maxje/llsoirl 

and E is the null space of G^ with E^E = I. 

Now we proceed. The methodology to solve the problem of (|28] ) is: 1) fix x, find the optimal P^; 2) 
fix Ps, find the optimal x. Based on this, we propose an algorithm to search for the optimal Pg and x 
as follows. 
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Algorithm 2: Take a feasible pj^^ as an initial point. Subsequently, find the optimal x^^) and then 
the optimal Ps^\ Then find the optimal x^^^ and so on. 

Remarks: We will see later that: for J = 2, the procedure always converges to the optimal Ps and x, 
while for J > 2 the procedure does not necessarily converge to the optimal solution. We will discuss 
this in the sequel. 

The algorithm |2] is not complete without providing the methods to find the optimal x for fixed Pg and 

find the optimal Ps for fixed x. Next, we provide such methods. 

First, we consider the problem: find the optimal Pg for fixed x. We need to solve 

a2 + P,|/ioP + (Po-i^s)xtR,,x 
max ■ — (31) 

a2 + max,-67{P,|ffojp + (^o - Ps)xtR^gx} 
s.t. G[Po°^-,Po]. 

Note that maxjg/{Ps|(7oj P + (Po — Ps)x^Rgx} denotes a polygonal line T whose vertices are located 
at the points Pg^k, A; = 0, 1, • • • , M, Pgfl = Pq™™, Ps,m = Po- An example is plotted in Fig. IH Note 
that the objective in (|3T]) is a linear fractional function and hence quasi-linear ll27l in each line segment 
Ps,kPs,k+i, A; = 0, 1, • • • , M — 1. Thus, the optimal Ps is one of the vertices of the polygonal, i.e., Ps^k, 
/c = 0, 1, • • • , M — 1. It is easy to find Ps^k- 




Second , we consider the problem: find the optimal x for fixed P^. We need to solve 

. a2 + P,|/io|2 + (Po-P,)xtR^x 

max min ■. — (32) 

X jei a2 + P,|5ojf + (Po-P.)xtR^^x 

s.t. x^x = 1. 
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We can show that the problem of (l32l ) is equivalent to 

max llvlP s.t. v^A,v < 1, 7 € / (33) 

where 

^3=K^'^^K^'\3^^ (34) 

R;, = (a2 + P,|/ioP)I + (Po-^s)Rh, (35) 

= [a" + P.|50,f )I + (^l - PsW,.3 e (36) 

Let the solution of ^ be v°. Then, the solution of ^ is given by x° = R^^/%7||R^ ^/%°||. The 
proof is given in Appendix |Hl 

The problem of (1331) belongs to the quadratic constrained quadratic programs (QCQP) which is in 
general difficult |[30l . |[3ll . Denote Z = vv^^ which enables us to rewrite it as 

maxTr(Z) s.t. Z ^ 0, rank(Z) = 1, Tr(AjZ) < 1, j G /. (37) 

Dropping the constraint rank(Z) = 1, the semi-definite relaxation (SDR) of the problem of (|37] ) is given 
by 

maxTr(Z) s.t. Z ^ 0, Tr(AjZ) < 1, j G /. (38) 

This is a semi-definite program (SDP) and can be effectively solved by CVX |32]. 

If the problems of (1331 ) and (1381 ) achieve the same objective value, we say the SDR of (1381 ) is tight 
(its solution Z° does not necessarily have rank one). Obviously, if the solution Z° of (1381 ) has rank one, 
the SDR (|38] ) is tight and the optimal solution v° of the problem of (|33] ) can be obtained by simply 
eigen-decomposing Z° = v°v°^. If Z° does not have rank one, the problem for the solution of (l33l ) is 
in general difficult |[30l . |[3T1 . But for the special case J = 2, the problem can be solved in polynomial 
time ||30l §2.2], lIST] Lemma 2.2.3]. We state it as the following theorem. 

Theorem 2: When J = 2, the SDR dJS]) is always tight and the optimal solution for the problem of 
f liil ) can be constructed in polynomial time. 

For J > 2, if Z° does not have rank one, we use Gaussian randomization procedure (GRP) to obtain an 
approximate solution based on the SDR solution |[28l . |[29l . In detail, we calculate the eigen-decomposition 
of Z° = UDU^ and generate v; = ^u/UD^/^^^ for Z = 1, • • • , L, where is a vector of zero-mean, 
unit-variance complex circularly symmetric uncorrelated Gaussian random variables, /U/ is chosen such 
that maxjg/ vj AjV; = 1 (i.e., the constraint of (l33l ) holds). Suppose 

/* = arg max ||vi|p. (39) 
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Then select v^'*) as an approximation solution. 

IV. Transmit Power Minimization under Secrecy Rate Constraint 
In this section, we address Problems |3] and ID 

A. DF-based protocol 
Problem [3] is to solve 

mill Pq = + ||w|p (40) 



Cj2 +P^|SrOp + wtRgW 







We provide a closed form solution to this problem which reveals that the optimal weight vector w° is a 
linear combination of h and g. We give the main result of the problem as a theorem. The proof is given 
in Appendix m 

Theorem 3: When 4^"|5oP - \ho? + 0, 

(CO) if C>^'o'"'",A2>-l 



Ai 

where C = (4^° — l)(T^/(|/ioP — 4'^"|5(op); Ai > 0, A2 < are the only two nonzero eigenvalues of 
R = (R/i — 4^°Rg)/(4^" |(/oP ~ l^op) yvith associated eigenvectors ui and U2 (see Lemma\I}. 
When 4^° I (/op ~ l^op = 0, the solution is 

po ^ pmin (.42) 



w 



(4«^'-l)aV6vi (43) 



where Ci > 0, ^2 < are the only two nonzero eigenvalues o/(R/i — 4^°Rg) with associated eigenvectors 
vi and V2 ( see Lemma |7]). 



Before ending of this subsection, we should point out that the results in |[22ll - ll24ll address only the 
case 4^=|5(op — |/iop > 0, thus our analysis here is more complete. Further, our analysis reveals that 
(R/i — 4^''Rg) has only two nonzero eigenvalues (one positive, and one negative) and, unlike ll22l - ll24l . 
provides simple expression for them. 
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B. CJ-based protocol 
Problem |4] is to solve 

mill Pq = P, + ||w|P (44) 

s.t. !— -J ^^-J = Pg > 0. 

wtR/jW + fj^ wtR^w + o"2 Ps ^ 

By denoting ||w|p = 7, x = w/||w||, vi = h/||h||, V2 = g/||g||, the problem of (l44l) can be rewritten 

as 

min Po = n + 7 (45) 
sT 2^|rf ^g5^.P.>0. 

7||hPxtviv{x + fj^ 7||gPxtv2V2X + fj^ -^s 

Let X be a feasible point, and xWivjx = z, z ^ [0, 1]. For fixed z and 7, a larger XW2V2X results 
in a smaller P^- With this and from Lemma |2j we know that the optimal XW2V2X equals F{z) where 
F{z) = 1 — (p-v/1 — z — — p'^)z)'^, p = |v|v2|. With this in mind, we can rewrite (|44l) as 

min Po = P, + 7 (46) 

1 1 1 

S.t 



^aiz + a2 '~iazF{z) + ai P<j ' 
[0,1], P, >0, 7>0 

where ai = ||hf (2^" - l)/\h^\\ 02 = ^t'(2^" - l)/\ho\\ as = ||gf (2«" - l)/(2^"|5on and = 
cr^(2^° — l)/(2^=|(7op). Further, we can rewrite the problem of (l46l) as 

. (7ai2; + Q2)[7a3i^(2;) + "4] , 
mm Po = — — — + 7 (47) 

7.2 ^[a^r [z) — aiz) + a4 — 02 

s.t. 7(a3P(z) — aiz) + — 02 > 0, 

ze [o,i],7>o. 

Before proceeding, we give a suboptimal solution which is turns out to be the same as the the suboptimal 
solution of ll24l . Please see Appendix |j] for details. 

Lemma 6: When z = is fixed, a suboptimal solution is pbtained as 

-fo.sub = -Ps.sub + 7sub (48) 
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where 

_j L if i;>l + V«3(l-/^2) 

I Q3(l-P=^) 

1/02 - l/[7suba3(l - /o^) + "4] ' ^^^^ 
Now we proceed. The methodology to solve the problem of (|471 ) is: 1) fix z, find the optimal 7; 2) 

fix 7, find the optimal z. Based on this, we propose an algorithm to search for the optimal 7 and z as 

follows. 

Algorithm 3: Take a feasible point z^^^ as initial point. Subsequently, find the optimal 7^^^ and then 
the optimal z^'^\ Then find the optimal 7*^^^, and so on. The procedure converges to the optimal 7° and 
z°. 

The algorithm [3] is not complete without providing the methods to find the optimal z for fixed 7 and 
find the optimal 7 for fixed z. Next, we provide such methods. 

First , we consider the problem: find the optimal 7 for fixed z. This corresponds to an optimization 
problem of a single variable 7, and the maximum is achieved at one of the following points: and the 
points with zero derivative. With this, we obtain 



7 = max<-!^ — ,0^ (51) 

y a^F{z) - aiz J 

where fi{z) = {a2a'j,F{z) — aia^^zY /[a'iF{z) — aiz + aia'izF{z)\. We can obtain 7 explicitly as a 
function of z (in this sense, we in fact reduce the original problem to a single variable optimization 
Po = Pf){z),z e [0, 1] explicitly). 

Second , we consider the problem: when 7 is fixed, find the optimal z. Let us denote the left side of 
the first constraint in (l46l ) by f^yiz), namely 

M^) = ■ (52) 

We know that the optimal z must maximize /^(z) which results in the minimal P, (see (|46b). The 
derivative of f^{z) given by 

/' (z)^^= _^ lasF'jz) ^^^^ 

dz (701^ + 02)^ {■jasF{z) + a^y 

It is easy to verify that /^(O) > 0, /^(l) < 0. Thus, z = is not the optimal point, and the optimal z 

must be the points with zero derivative. Note that Ps > (i.e., f^{z) > 0) holds only when 701^ + 02 < 

jasF{z) + 04 which determines an interval {z_,z) C [0,1]. Here z and z can be expressed in closed 

form from the fact: if F{z) = f3iz + /32 has real roots over [0, 1], then its roots can be expressed as 
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z = 1/(1 + ul) Where uq = (pTT^i ^p^{l - p^) + (/32 - 1 + /92)(p2 _ _ /32))/(/32 - 1 + /o^). 
With these, we can restrict our attention to the root of ( [53] ) over (z, z). 

To proceed, we need the following result. The proof is given in Appendix iKl 

Property 2: For f^{z) defined in (l52l ). ^^|^, < 0/or the stationary point z' G {z_,z) (i.e., the point 

According to Property |2j we know that the equation (1531 ) has a unique root z' such that when z < z' , 
fif{z) > and when z > z', f!^{z) < 0. This property ensures that the Newton method would be very 
effective in searching for z' and would enjoy quadratic convergence. 

V. Numerical Simulations 

In this section we provide some numerical simulations to illustrate the proposed solutions. We use the 
same system configuration as that in |[24l . where source, relays, destination and eavesdroppers are placed 
along a line. Channels between any two nodes are modeled as a line-of-sight (LOS) channel pQd~'^^'^e^^ , 
where d is the distance between two nodes, po is a constant, c is the path loss exponent, and 6 is the 
phase uniformly distributed within [0, 2tt). In our simulations we set c = 3.5 and pQ = 1. We assume the 
distances between relays are much smaller than the distances between relays and source/desination, such 
that the path loss between different relays and source/destination can be taken as approximately the same. 
Similarly, the path loss between different eavesdroppers and source/destination/relay are approximately 
the same as well. The results are obtained using Monte-Carlo simulations consisting of 500 independent 
trials. 

First, we vary the position of the destination so that the source-destination distance changes from 10 m 
to 100 m, as shown in the upper row of Fig. |3] The source-relay distances are fixed at 5 m, the number 
of relays is = 10, the source-eavesdropper distances are fixed at 50 m, the power constraint is fixed at 
30 dBm, the secrecy rate constraint is fixed at 1 bits/s/Hz. The secrecy rate for a single eavesdropper and 
multiple eavesdroppers is depicted in Figs. |4] and [51 respectively. From these two figures, one can see that 
when the destination moves past the eavesdropper direct transmission cannot sustain positive secrecy rate. 
On the other hand, both DF and CJ maintain positive secrecy rate even when the destination is further 
away from the source than the eavesdropper. The fact that there is a cooperation advantage even when 
the destination is at the same location as the eavesdropper is because of the phases differences of the 
corresponding channels. Although the propagation environment would be the same for both destination 
and eavesdropper in that case, the phases will be different due to different receiver phase offsets. The DF 
scheme yields the higher secrecy rate, while the optimal and suboptimal CJ schemes produce the same 
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average rate. 

Similar observations can be drawn from Fig. [5] for the case of multiple eavesdropper as far as the 
advantage of cooperation over direct transmission is concerned. 

In Fig. |6l the secrecy rate for a fixed configuration and variable number of eavesdroppers is shown. It 
can be seen from Fig. [6] that, when the number of eavesdroppers increases, the suboptimal solution for 
DF in Lemma [5] becomes inferior as compared to the optimal solution. The minimal transmit power is 
depicted in Fig. |7] For comparison purposes, the suboptimal solution for CJ in Lemma |4] and the direct 
transmission result are also shown on the same figure. 

Second, in Fig. [H we show the performance when the eavesdroppers' positions change while the 
source-destination distance is fixed at 50 m and the source-relay distances are fixed at 5 m. When the 
source-eavesdropper distance changes from 25 m to 100 m as shown in the lower row of Fig. [3j the 
minimum transmit power for CJ first increases a little, and then decreases, while the minimum transmit 
power for DF always decreases. The results show that cooperation can significantly improve the system 
performance as compared to direct transmission. In particular, when the source-eavesdropper distance is 
smaller than 65 m, using direct transmission there is no level of transmit power than can meet the secrecy 
rate constraint. Also, for source-eavesdropper distance greater than 85 m direct transmission and the CJ 
scheme are equivalent in terms of the minimum required transmit power. The DF approach requires 
significantly smaller power to meet the secrecy rate constraint. It is interesting to note that in the average 
sense, the suboptimal solution for CJ in Lemma |4] and Lemma [6] is a very good approximation of the 
optimal solution. 

VI. Conclusion 

We have given explicit constructions for the optimal relay weights and source transmission power for 
maximizing the secrecy rate or minimizing the total transmit power (source transmission power plus 
relay power) under secrecy rate constraint using the DF and CJ protocols in the presence of a single 
eavesdropper or multiple eavesdroppers. We present numerical results to compare the secrecy rate under 
our optimal solutions with the secrecy rate under the sub-optimal solutions in |[22l - ll24ll . Numerical 
results illustrate that cooperation can significantly improve the system performance as compared to direct 
transmission. 
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Appendix A 
Proof of Lemma [U 

Let A 7^ be the eigenvalue of rr^ — ss^ associated with the eigenvector a. Thus we have (rr^ — sst)a = 
Aa which leads to a = (1/A)[r(r1^a) — s(s1^a)]. Thus, a has the form of a linear combination of r and s. 
With this, we let a = vrir + 7r2S where tti and 1^2 will be determined as follows. Since a has unit norm, 
we have 

|7rip||r|p + 7ri7r2(s'''r) + 7ri7r2(r^s) + |7r2p||s|p = 1. (54) 
On the other hand, by inserting a = vrir + 7r2S into (rr^ — ss1^)a = Aa, we get 

[vridlrf - A) + ^2(rts)]r = [^i(str) + 7r2(||sf + A)]s. (55) 
Since r and s are linearly uncorrected, we have 

TTidlrf - A) + ^2(r^s) = 7ri(str) + ^2(||sf + A) = (56) 

which leads to 

A2 _ (||r||2 _ ||s||2)A - (||rf ||sf - |rts|2) = (57) 
A - llrlP 

712 = 1 TTi. (58) 

rTs 

From Cauchy inequality, we get ||r|p||s|p — jr^'sp > 0. Thus, The equation (1571 ) has a positive root and 
a negative root. As a result, we obtain A and the corresponding vri and 1^2- 

Appendix B 
Proof of Lemma [2] 

The solution z° of (fT4l) is a linear combination of di and d2, which follows from its optimality 
condition d2d2Z — ^ididjz — /i2Z = or further /i2Z = (d2z)d2 — (/iidjz)di where /ii and /i2 are 
Lagrange multipliers. Note that e^'^^z" = e~'^^cidi + |c2|d2 is also solution of (fT4b . where 6*2 is the 
argument of C2. Consequently, we can restrict C2 > 0. Inserting z = cidi + C2d2 into the constraints and 
objective, results in 

|cip + c| + c^C2d|d2 + ciC2d5di = 1 (59) 
|ci|^ + c^|d|d2p + CiC2d|d2 + ciC2d5di = q (60) 
zMad^z = 1 - |ci|2(l - |d|d2|2). (61) 
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From (|6T]) . we need to minimize |cip. From ( [59l ) and (l60l) . we get c^l — |d|d2p) = 1 — Q which leads 



to C2 = ^/(l — — '"^)- By denoting ci = |ci|e'^ where 9 is the argument of ci, we can rewrite (l59 
as 



|ci|2 + |ci|c2r(e-(<^+^) + e^('^+^)) + (c^ - 1) = 0. 



It is not difficult to show that the optimal 9 given by 



vr ■ 



if ci - 1< 
if c| - 1 > 



and the optimal |ci| is given by 



ci 



y/q- C2r if C2 - 1 < 
^2?^ - a/^ if c| - 1 > 
With these, we obtain the optimal ci = (c2r — ^)e^'^'^~'^\ Further, from (|6TI ). we obtain 

zM24z = 1 - (C2r - - = 1 - (ryr^ - ^(1 " r^)^)^. 



(62) 



(63) 



(64) 



(65) 



Appendix C 
Proof of Theorem [T] 



First, we derive the optimal weight vector w for fixed Pg > P™^"^. By denoting = y/Po — Pg x, we 
can rewrite the problem of (fTSl ) as 



1 , (a^ + Ps\ho\^ + (Po - Ps)||h||2xtu24x' 
max - log ' 



+ PsboP + (Po - P.)||gPxtuiulx 



(66) 



s.t. x^x 



1. 



Let X be a feasible point and x^uiujx = z, z £ [0, 1]. For fixed z, a larger x^U2U2X results in a larger 
objective value in the problem of ( [66l ). With this and from Lemma |2j we know that x^U2U2X equals 



L{z). Thus, the problem of (l66l ) can be rewritten as 



, 1 1 f^^ + PsM'' + (Po - Ps)||hf 
max M(2;) = - log ' 



+ Ps\go\^ + {Po - Ps 



(67) 



s.t. < z < 1. 



This is an optimization problem of a single variable z. It is easy to show that M'(0) > 0, M'(l) < 0. 
Thus, the optimal z must be the points with zero derivative, i.e., M'(z) = 0. As a result, the solution of 
(|67] ). denoted by z(Ps) as a function of Pg, can be expressed in closed form (l20l ). 
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Next we consider the optimal Pg and let it be P°. We can state that P° is also the solution of the 
following associated problem 



To see why this is the case, let assume the solution of (1681 ) is P^ but not P°. Denote the objective in (|67] ) 
as ^log{Ji{z,Ps)). Note that the objective in (l68]l is exactly i log( Ji(z(P°), P,)). Recall that for fixed 
P^, the solution of ^ is which leads to i log(Ji(z(P,'), -Pi)) > l^og{Ji{z{P°), P;,)). On the 

other hand, ^log{Ji{z{P°), P°)) > ^ log(Ji(z(Pi), P,')) since P° is the solution of ((B]). Combining 
both gives ^ log(Ji(2;(P°), P°)) > ^ log(Ji(2;(P°), P^')) which violates the assumption that the solution 
of dM) is Ps but not P°. 

Further, Ji{z,Ps) is a linear fractional function known to be quasi-linear jlTl . thus the maximum 
always occurs at one of the two ends. Thus, the solution of ([T5] ) must be Pq"™ or Pq. With this, we can 
also obtain the optimal w according to Lemma [2] 



From the objective in the problem of ( |23] ). it should hold that (Pq— Ps)aiz+a2 < {Po — Ps)ct3G{z)+a4 
for some Ps G [0,Po], z G [0, 1]. In other words, it should hold (Pq — Ps){a3G(z) — aiz) > Q2 — 04 
for some Pg G [0,Po], z G [0,1]. Denote K{z) = a3G{z) — a\z, z G [0,1]. It is easy to verify that: 
1) ir(0) = 03(1 - r/2) > 0; 2) K"lyz) < 0; 3) K'{z) ^ +00 as z ^ 0, K'{1) < 0. Here K'{z) and 
K"{z) denote the first- and second- order derivatives, respectively. With these, first, if 02 < 04, then 
(Pq — Ps)K{z) > 02 — 04 for some Ps G [0, Pq], z G [0, 1] holds since i^(0) > 0; second, if 02 > 04, 
we know that K{z), z G [0, 1] achieves its maximum at zq with K'{zq) = (i.e., the unique root of the 
equation a^G'{z) — oi = 0), and (Pq — Ps)K{z), Ps G [0, Po],z G [0, 1] achieves its maximum PoK{zo), 
then, the condition should be PoK{zo) > 02 — Q!4. 




(68) 



s.t. P, G[Po-'",Po]. 



Appendix D 



Proof of Lemma [3] 



Appendix E 



Proof of Lemma |4] 



When z = is fixed, the problem of (1231 ) is reduced to 




(69) 



s.t. P, G[0,Po]. 
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This is an optimization problem of a single variable Pg, and the maximum is achieved at one of the 
following points: 0, Pq and the points with zero derivative. After some calculations, the desired results 
are obtained. 

Appendix F 
Proof of Property [T] 

According to (1231 ). we can write 

i?,(z) = log(l + ^) - log (1 + \ ) (70) 

qiz + q2 q-iG{z)+qi' 

where qi = {Pq - Ps)ai/ Pg, q2 = 012/ Ps, qz = (-Po - Ps)a3/Ps and qi = a4/Ps. It follows from 

dRs I ^ qi 

<9z {qiz' + q2 + l)iqiz' + q2) 

+ = (71) 

^ {q,G{z') + q, + l){qsG{z') + q,) ^ ' 



that 

iq^Giz') + q^ + l){q3G{z') + q^) q^G' {z') 



{qiz' + q2 + l){qiz' + q2) qi 

On the other hand, we know 

d^Rsi qli2{qiz' + q2) + 1) 



(72) 



9^2 I.' l^q^z' + q2 + l)iqiz' + q2)]'^ 

iq3G'iz'))H2iq,Giz')+q^) + l) 
[(q,Giz')+q^ + l)iq3G{z')+q^)]^ 

q3G"{z') 

+ {q^G{z') + q^ + l){q3G{z') + q4)- 

Inserting ( 1721 ) into ( 1731 ) and using the fact: G"{z') < and q^G{z') + ^4 > giz' + q2 (since z' G (z, z)) 

leads to the desired result. 



Appendix G 
Proof of Lemma [5] 

When the constraint w^G = is added, the problem of (l27l) is reduced to 

1 a2 + P,|/io|2 + wtR;,w 

o^og 2 I p ri 12T 

p,,w 2 + Psm.a-Kj(zj{\goj\^\ 



s.t. P, G [Po™", i^o], wtw = Po - P„ wtQ = 0. 
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From w^w = Pq — Pg and w^'G = 0, we can obtain w = ^/Pq — Ps Ez and z^^z = 1 which, when 
inserted into (1741) . results in 

1 , + Ps\ho\^ + (Po - Ps)ztEtR;,Ez 
max - log 5— — -J- / :>) 

P,,z 2 + Psma.yij(zi{\goj\^} 

s.t. G [Po'"'^Po], ztz = l. 
By using the fact z^E^^R/jEz achieves its maximum at z = E^h/||El^h||, we rewrite the problem of (175] ) 



as 

1, ^^ + P.|/^oP + (Po-Ps)||Ethf 

max - loff 7^ TTTT (76) 

ft 2 ^ (T2 + P,max,eHl50jf } 

s.t. P,G [Po'^^Po]. 

The objective in (1761 ) is ^ log /2(Ps). Note that /2(Ps) is quasi-linear ll27l . thus the maximum occurs at 
one of the two ends, namely, p™™ or Pq. 

Appendix H 
Proof of Equivalent Problem (i33l) 

From the constraint x'I'x = 1, we can rewrite + Ps|/ioP = (d^ + Ps\ho\'^):x) 'k and cr^ + P^j^ojP = 
((T^ + Ps|(7oj p)x'l"x which enables us to rewrite the problem of (l32l) as 

x'l'R/jX 

max min — — — . (77) 
l|x||=i jei xtR^x 

Further, by denoting u = R^^^x/||R^^^x||, we can rewrite the problem of (1771) as 

max min —j— s.t. u^u = 1. (78) 

u jei uTAjU 

By introducing the slack variable y to rewrite min^g/ l/(u^Aju) = 1/y, we can rewrite the problem of 
(Ell) as 

max — s.t. u^u = 1, u^AjU <y,j G /. (79) 
y y 

Let V = u/y^, then u^u = 1 is equivalent to ||v|p = 1/y and we can rewrite the problem of (1791 ) as the 
problem of (l33l ). Let the solution of (l33l ) be v°, then the solution of (1781 ) is u° = v°/||v° ||, the solution 
of dTTI) is x° = R^^/^v7||R^^/%°||. 
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Appendix I 
Proof of Theorem [3] 

The first constraint of (l40b leads to 

(4«°l5oP - \ho?)Ps = wt(R, - 4^=R,)w - (4«= - l)a\ (80) 

There are two cases. If 4^=|g'oP - \ho\^ / 0, it follows from that Pg = w^^Rw + C- With this, we 
can rewrite the problem of (l40l ) as 

min wW + w^Rw + C s.t. w^Rw + C > -Pn™''. (81) 

w 

The optimal relay weight vector w° has the form of (diUi + (i2U2) which follows from its optimality 
condition of (ISTT) ETl 

w + Rw - z/Rw = 0. (82) 

where i/ > is the Lagrange multiplier. Indeed, inserting R = Aiuiuj + A2U2U2 into (l82l ) gives 
\v = {v — l)Ai(uJw)ui + {v — l)A2(u|w)u2. With this, we can rewrite dST] ) as 

min (1 + \i)\di\^ + (1 + \2)\d2\^ + C (83) 

di,d2 

S.t. Ai|di|2 + A2|d2|^ + C> JT""- 

There are several case. If C > ^T'"' ^2 > -1, then the solution is |fii| = 0, jdal = 0; If C > iT'""' 
A2 < -1, then the solution is \di\ = 0, |(i2p = (C - ^o'''°)/|A2|; If C < ^o""'", the solution is \di\^ = 
(Po™--C)/Ai, 1^2! =0. 

Similarly, we can obtain the other results in Theorem |3] 

Appendix J 
Proof of Lemma [6] 

When z = is fixed, the problem of ( |47] ) is reduced to 

„ Oi2hc(3(l - p'^) + a^] ^^^^ 

mm Po = V, 2^ + ^ (^4) 

7 7a3(l — p^j + 04 — a2 

Q^2 ~ CK4 

03(1 - p^) 

This is an optimization problem of a single variable 7, and the maximal is achieved at one of the following 
points: and the points with zero derivative. After some calculations, the desired results are obtained. 
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It follows from 



Appendix K 
Proof of Property [2] 



^^''^ - "^"^ ^tf!"'^ .. =0 (85) 



that 



dz (7012;' + 02)^ (7Q3F(z') + 04)2 

(7Qiz' + a2)^ ai 

On the other hand, we know 

92/7. I 2(7ai)2 2{^a,F'{z')Y 



(86) 



^-,2 \z' {^aiz' + {lOi-iF{z') + a^Y 

jasF"{z') 

^ {jasF{z') + a,r- ^ ' 

Inserting (l86l) into ([87]) and using the fact: F"{z') < and 703^(2') +04 > ^aiz' +a2 (since z' G (z, z)) 
leads to the desired result. 
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Fig. 3. Simulation model: source, relays, destination, eavesdroppers are placed along a Une; Upper one for Fig. 4, 5, 7; Lower 
one for Fig. 8. 




Fig. 4. Secrecy rate vs. source-destination distance (power constraint: 30 dBm, source-relay distance: 5 m, number of relays: 
iV = 10, one eavesdropper (J = 1), source-eavesdropper distance: 50 m). 
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Fig. 5. Secrecy rate vs. source-destination distance (power constraint: 30dBm, source-relay distance: 5 m, number of relays: 
A^^ = 10, number of eavesdroppers: J = 7, source-eavesdropper distance: 50 m). 
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Fig. 6. Secrecy rate vs. number of eavesdroppers (power constraint: 30 dBm, source-relay distance: 5 m, number of relays: 
N = 10, source-destination distance: 25 m, source-eavesdropper distance: 50 m). 
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Source-destination distance (m) 



Fig. 7. Transmit power vs. source-destination distance (secrecy rate constraint: Ibits/s/Hz, source-relay distance: 5 m, the 
number of relays: A'^ = 10, one eavesdropper (J = 1), source-eavesdropper distance: 50 m). 
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Fig. 8. Transmit power vs. source-eavesdropper distance (secrecy rate constraint: 1 bits/s/Hz, source-relay distance: 5 m, 
number of relays: N = 10, one eavesdropper (J = 1), source-destination distance: 50 m). 
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